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ABSTRACT 

We present broadband observations and spectral modeling of PKS B0008-421, and identify it as an extreme 
gigahertz-peaked spectrum (GPS) source. PKS B0008-421 is characterized by the steepest known spectral 
slope below the turnover, close to the theoretical limit of synchrotron self-absorption, and the smallest known 
spectral width of any GPS source. Spectral coverage of the source spans from 0.118 to 22 GHz, which includes 
data from the Murchison Widefield Array and the wide bandpass receivers on the Australia Telescope Compact 
Array. We have implemented a Bayesian inference model fitting routine to fit the data with internal free-free 
absorption, single and double component free-free absorption in an external homogeneous medium, free-free 
absorption in an external inhomogeneous medium, or single and double component synchrotron self-absorption 
models, all with and without a high-frequency exponential break. We find that without the inclusion of a high- 
frequency break these models can not accurately fit the data, with significant deviations above and below 
the peak in the radio spectrum. The addition of a high-frequency break provides acceptable spectral fits for 
the inhomogeneous free-free absorption and double-component synchrotron self-absorption models, with the 
inhomogeneous free-free absorption model statistically favored. The requirement of a high-frequency spectral 
break implies that the source has ceased injecting fresh particles. Additional support for the inhomogeneous 
free-free absorption model as being responsible for the turnover in the spectrum is given by the consistency 
between the physical parameters derived from the model fit and the implications of the exponential spectral 
break, such as the necessity of the source being surrounded by a dense ambient medium to maintain the peak 
frequency near the gigahertz region. This implies that PKS B0008-421 should display an internal HI column 
density greater than 10 io cm 2 . The discovery of PKS B0008-421 suggests that the next generation of low 
radio frequency surveys could reveal a large population of GPS sources that have ceased activity, and that a 
portion of the ultra-steep spectrum source population could be composed of these GPS sources in a relic phase. 
Keywords: galaxies: active — galaxies: individual (PKS B0008-421) — radiation mechanisms: general — 
radio continuum: general — radio sources: spectra 


1. INTRODUCTION 

Gigahertz-peaked spectrum (GPS) radio sources are a class 
of active galactic nuclei (AGN) that have played a pivotal 
role in shaping our understanding of the environmental prop¬ 


erties and evolutionary paths of radio galaxies. These pow¬ 
erful radio sources are defined by a concave spectrum that 
peaks at gigahertz frequencies, steep spectral slopes on ei¬ 
ther side of the turnover, small linear sizes (~0.1 — lkpc), 
minimal variability over timescales varying from hours to 
decades, and low radio polarization fractions (lO’Dea et all 
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1199 ll) . Two closely related groups of radio sources are high 
frequency peakers (HFPs) and compact steep spectrum (CSS) 
sources, which display many of the same physical prop¬ 
erties as the GPS population but with different peak fre¬ 
quencies and linear sizes. HFPs peak above a gigahertz 
and have smaller linear sizes than GPS sources (5.0.1 kpc ) 

dDallacasa et al.ll2000fc iTinti et al .120051 lHancock et al.l2010h . 

while CSS sources peak below gigahertz frequencies and 
have larger linea r sizes than GPS sources (~ 1 — 10 kpc) 
dFanti et alJll990l dp). The hosts of HFP, GPS and CSS 
sources can be associat ed with quasar s, radio g alaxies, or 
Seyfert galaxies (jGelderman & Whittlelll994t Ide Vries et alJ 
Il997t IStanghellinill2003l) . A comprehensive review of GPS 
and CSS sources has been presented bv lO’Deal (119981) . with 
additional and more rec ent research summ aries provided by 
IStawarz et alJ (120081) and lMarr et al.l (120141) . 

Very long baseline interferometry (VLBI) imaging of HFP, 
GPS and CSS sources has allowed these radio sources to be 
further split into two distinct classes based on their morphol¬ 
ogy. The first morphological class is defined by a core-jet 
structure, and su ch sources are generally associated with hig h 
redshift quasars (tStanghellini et alii 19971 lOrienti et aPl2006l) . 
The second morphological class of GPS and CSS sources 
is characterized by dominant small scale structures that are 
reminiscent of the large scale radio lobes of powerful radio 
galaxies, with inverted or flat spectrum cores surrounded by 
two steep-spectra lobes. GPS and CSS sources that display 
these characteristics are referred to as compact symmetric ob¬ 
jects (CSOs), and are generally associated with low redshif t 
radio galaxies (IWilkinson et al'lfl994t Ide Vries et al.ll2009bl) . 
CSOs have been the focus of extensive study because such 
morphologies suggest that these s ources co uld be the young 
precursors to large radio g alaxies d Marr et al.ll2014l). _ 

Thi s ‘youth’ hypoth e sis dPhillins & Mutelll982HFanti et al.l 
119951: ISnellen et alJ 120001) has observational support 
from hotspot measurements (lOw sianik & Conwavl 119981 
iPolatidis & Conwavl 120031 iGugliucci et al.l 120051) and from 
high-frequency spe ctral break modeling dMurgia et al.ll 19991 
lOrienti et al.l 120101) . However, such an interpretation is still 
contentious, as statistical studies of the luminosity functions 
have demonstrated an over abundance of CSOs relative to 
the number of large r adio galaxies (lO’Dea & Bauml 119971 
iReadhead et al.l 119961 I An & Baaiil 120121) . An alternative 
hypothesis is the ‘frustration’ model, which implies that 
these sources are not young but are confined to small 
spatial scales due to unu sually high nuclear plasma density 
(Ivan Breugel et al.l 11984 ). It is possible that both these 
scenarios may apply to GPS sources, as young sources with 
constant AGN activity could break through a dense medium 
given sufficient time, whil e sou rces with intermittent AGN 
activity may stagnate (An & Baanll2012l) . 

One of the reasons that there has not been a resolution be¬ 
tween these two competing hypotheses is because the absorp¬ 
tion mechanism responsible for the spectral turnover still re¬ 
mains an open area of debate. Synchrotron self-absorption 
(SSA) by the relativistic electrons internal to the emitting 
source, and free-free absorption (FFA) via an external ionized 
screen of dense plasma, are two comm only proposed mod¬ 
els for the turnover dKellermannl 1966 ). Both these absorp¬ 
tion models replicate some physical features of GPS sources 
but fail to explain others. For example, the observed correla¬ 
tion between the turnover frequencies and linear sizes of GPS 
sources is well justified within the SSA framework, while 
FFA via a homogeneous medium can not replicate such a re¬ 


lationship (IO’Dealll998l) . However, FFA models that invoke 
an external inhomogeneous medium where electron density 
decreases with distance f rom the radio jet (e.g. lBicknell et all 
IT997l:lKuncic~et~aT]IT998i) . or dense ionized interstellarclouds 
that co-exist with the rela t ivistic electrons (e.g . iBe gelmanl 
119991 IStawarz et alJ 120081 iMaciel & Alexanden 120141) . are 
able to recreate this correlation. While it has been shown 
that it is likely that this ‘frustration’ sceanrio is inconsis- 
tent with the properties of many GPS sources (ISnellen et alJ 
l2000Ll2003l:lde Vries et al.l2009al) . there is mounting evidence 
from many observational studies of individual GPS sources, 
that have morphologies consistent with CSOs, demonstrating 


Kameno et all 2000; 

Marr et alJ 2001; Orienti & Dallacasa 

2008; Tremblay et al. 

2008]; Marr et ah. 20141; ITingav et al. 


2015). 


It is possible to differentiate between FFA and SSA models 
through comprehensive statistical fitting of GPS spectra, pro¬ 
vided the spectra are well s a mpled below the turnover. For 
example, iTingav & de Kooll (120031) performed a detailed in¬ 
vestigation of different absorption scenarios via model fitting 
to the radio spectrum of the GPS source PKS B1718-649. 
They concluded that SSA was the most likely contributor to 
the inverted spectrum, but emphasized that thejnodified in¬ 
homogeneous FFA model of iBicknell et al] (119971) could not 
be excluded. The reason the degeneracy between SSA and 
FFA models has remained unresolved is because these past 
investigations lacked comprehensive spectral coverage below 
the turnover, where the distinction between the different ab¬ 
sorption models becomes pronounced. With new low radio 
frequency t elescopes such as the Murchison Widefi eld Ar¬ 
ray (MWA; iLonsdale et al.ll2009l : ITingav et al.ll2013l) a nd the 
LOw-Frequency ARray (LOFAR; Ivan Haarlem et al.l 120131) 
becoming operational, astronomers now have unprecedented 
frequency coverage below the turnover, allowing for a more 
comprehensive spectral comparison of the different absorp¬ 
tion models. Additionally, the new wide bandpass receivers 
installed on existing radio telescope s, su ch as the Com pact 
Array Broadband Backend (CABB; IWilson et al.l 1201 ll) on 
the Australia Telescope Compact Array (ATCA), provide new 
broadband spectral coverage at gigahertz frequencies. Such 
broadband spectral coverage is also vital to understanding 
the nature of high-frequency spectral breaks in GPS sources, 
which can heavily influence spectral fits below the turnover 
(ITingav et al.ll2015l) . 

In light of these technological advancements, we present 
MWA and CABB observations of GPS source PKS B0008- 
421 (RA = 0 0:10:5 2.5, Dec. = -41:53:10.6 (J2000); 2 = 
1.12; lLabiano et all l200l . PKS B0008-421 was selected 
for this study because it has the steepest spectral slope be¬ 
low the turnover for sources obs erved during the MWA Com - 
missioning Survey (MWACS; iHurlev-Walker et ahl l2014t) . 
P KS B0008-42 1 was identified as a candidate GPS source 
by lO’Dea et al.l (1 9911) . While several studies have demon¬ 
strated properties of PKS B 0008-421 that are consistent with 
a GPS source classification (lLabiano et al. 2007; Jacobs et alJ 
120TI1) . we present the first flux density measurements below 
the turnover that confirm this characterization. VLBI mea¬ 
surements of PKS B0 008-421 reveal a morphology t hat is 
consistent with a CSO (lKinglll994l : lJauncev et al.ll2003l) . Fur¬ 
thermore, PKS B0008-421 has extensive temporal informa¬ 
tion and wide spectral coverage from 0.118 to 22 GHz, since it 
has acted as a primary or secondary calibrator for most south¬ 
ern hemisphere radio telescopes, allowing us to overcome the 
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epochal and spectral limitations of previous studies of GPS 
sources. 

The purpose of this paper is to present unprecedented spec¬ 
tral coverage of this extreme GPS source, and to demonstrate 
that advances in wide bandpass receivers and low radio fre¬ 
quency coverage now allow us to place stringent observational 
constraints on the absorption models and physical environ¬ 
ments of GPS sources. The relevant data reduction proce¬ 
dures performed for the MWA, CABB and archival observa¬ 
tions are discussed in §[2] In §0 the Bayesian fitting routine 
that we implemented to assess the different absorption model 
fits to PKS B0008-421 is outlined. Relevant features of the 
absorption models, and their respective fitting statistics, are 
presented in §Q] In §0the impact the absorption models fits 
have on our understanding of the environment of PKS B0008- 
421, and the absorption models used to describe GPS sources 
as a whole, are discussed. 

2. OBSERVATIONS, DATA REDUCTION AND ARCHIVAL DATA 

PKS B0008-421 was regularly observed between 1969 and 
2014 by a variety of radio telescopes. Most of the observa¬ 
tions have been performed by the Parkes 64 m radio telescope, 
the ATCA, the Molonglo Observatory Synthesis Telescope 
(MOST) and the Giant Metrewave Radio Telescope (GMRT). 
Combining all these observations with the MWA and ATCA 
observations provides a spectral coverage of PKS B0008-421 
from 0.118 to 22 GHz. These observations are summarized in 
Tableland plotted in FigureQ] PKS B0008-421 was spatially 
unresolved in all these observations. 

2.1. ATCA Observations and Data Reduction 

As part of project C007 (PI Stevens), PKS B0008-421 
was observed by the ATCA in the 6C array configuration on 
2010 June 9-10. The observation was conducted using the 
CABB backend system, which gives two instantaneous 2 GHz 
bandwidth for both linear polarizations, at central frequencies 
5.5 GHz and 9.0GHz for a total integration time of 32 min¬ 
utes. Each observation used 1 MHz channels and a 10 sec¬ 
ond correlator integration time. PKS B0008-421 was also ob¬ 
served for project C2697 (PI Ryder) with the ATCA in the 
6D configuration on 2012 June 14-15. This observation was 
conducted at a central frequency of 2.1 GHz and for a total 
integration time of 121 minutes, using the same bandwidth, 
channel size and correlator integration time as the 2010 obser¬ 
vation. PKS B1934-638 was used for flux density, bandpass 
and phase calibration for both observations. 

The data for both obser vatio ns were reduced using the 
MIRIAD software package (lSault et al.lll995 ). The known re¬ 
gions of radio frequency interference (RFI) and lower sensi¬ 
tivity in the CABB system were initially flagged. RFI was 
excised for PKS B1934-638 using the automatic flagging op¬ 
tion in pgflag and manually with biflag. The bandpass, 
gain and leakage solutions were estimated using PKS B1934- 
638. The gain calibration was performed over four 512 MHz 
wide frequency bins and ten second time intervals. The cali¬ 
bration solutions were transferred to PKS B0008-421 and the 
flux density scale was set to that defined by PKS B1934-638. 

The flux density of PKS B0008-421 was measured in every 
spectral channel of the visibilities using the task uvfmeas. 
The uncertainty of each flux density measurement is domi¬ 
nated by thermal noise and gain calibration errors. The ther¬ 
mal noise was estimated by dividing the visibility space into a 
large number of independent cells of a size given by a window 
function, and by using the physical antenna size, bandwidth. 


total system temperature, time averaged number of baselines 
in the cell, and total observation time (Eqn. 11; iMoralesl 
120051) . The therma l noise is Gaussian and independent be - 
tween the channels dBowman et alJl2006l : IWilson et alJl201 lh . 
The gain calibration errors were estimated by examining the 
variation of flux density of the source and calibrator over dif¬ 
ferent baselines and by comparison of the root-mean-square 
residuals between the measured visibilities and the model vis¬ 
ibilities of a point source. The gain calibration errors are inde¬ 
pendent since the CABB continu um c hannels have a ‘square’ 
response (Fig. 12: iWilson et al.l 120 1 11) but may not be nor¬ 
mally distributed since they are antenna dependent. However, 
the magnitude of the deviation from a Gaussian distribution 
is dependent on the uniformity of the gain response of the an¬ 
tennas in the array dBoonstra & van der Veenl2003i) . Since the 
gain response of the antennas of the ATCA are within ~ 1% 
of each other (IWilson et al.H201 II) . it can be assumed that the 
uncertainties are approximately normally distributed. Addi¬ 
tionally, the effect of confusing sources on the spectrum of 
PKS B0008-421 is negligible, with the next brightest source 
within the full width at half maximum of the primary beam 
of flux density between ~ 0.4 to 0.6 mJy for all observing 
frequencies. For the same frequency range, PKS B0008-421 
varies between ~ 4 to 0.5 Jy. 

The flux density measurements and their associated uncer¬ 
tainties were grouped into bins of approximately 200 MHz, 
with variations in bin size occurring since sections of the 
spectrum had been flagged due to RFI. Images were also pro¬ 
duced at 1.8, 5.5 and 9.0 GHz to ensure that the visibilities of 
PKS B0008-421 were not contaminated by confusing sources 
and to confirm that PKS B0008-421 was spatially unresolved 
at all frequencies. 

2.2. MWA Observations and Data Reduction 

The MWA observed PKS B0008-421 during MWACS, a 
survey which covers approximately 21 h < RA < 8 h, — 55° 
< Dec. < — 10° over t hree frequency bands center ed at 0.199, 
0.150 and 0.180GHz (iHurlev-Walker et al.ll20l4h . MWACS 
was conducted between September and December 2012 by 
taking drift scans at declinations approximately — 27° and 
— 47°, using two different 32-tile subarrays of the full 128- 
tile array (where a tile consists of 16 dual-polarization dipole 
antennas). The survey has a 3 arcminute resolution and noise 
level of ~ 40mJy beam -1 at 0.150 GHz. 

Given that the MWA flux density measurements are vi¬ 
tal in constraining the spectral slope below the turnover of 
PKS B0008-421, it is important to summarize the absolute 
flux density calibration that was undertaken for MWACS. The 
absolute flux density scale was empirically evaluated by boot¬ 
strapping from a source that had been observed by northern 
and southern hemisphere radio telescopes. The source 3C32 
was selected since it has wide spectral coverage that is well 
fit by a simple power-law, shows no evidence of variabil¬ 
ity, and is unresolved in MWACS and in the source catalogs 
used to fill in its spectral energy distribution. The flux den¬ 
sities of sources in the Dec. — 27° maps were corrected by 
the least-squares fit to the catalog flux densities of 3C32. The 
Dec. — 47° mosaic was corrected to the flux density scale of 
the Dec. —27° mosaic by using «600 unresolved sources 
in the overlapping region of the two mosaics (see §3.4 of 
Hurley-Walker et al. 2014 for more details). PKS B0008- 
421 was unresolved in all images that contributed to MWACS. 
Note that we have used all three flux dens ities measur ed by 
MWACS, while the published catalog of lHurlev-Walker et al.1 
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Table 1 

A summ ary of the observations of PKS B0008-421 used in the spectral modeling. Note that the reported tWillsl 119751 fluxes have been corrected to the 
IBaars et all 119771) flux density scale. The epochs of the observations are presented as accurately as possible, but when only a broad timeframe is known, it is 

presented in the format YYYY-YYYY. 


V 

(GHz) 

Hu 

(Jy) 

<J 

(Jy) 

Epoch 

Telescope 

Reference (Survey) 

0.118 

0.7 

0.2 

2012 Dec 01 

MWA 

Hurlev-Walker et al. (2014) (MWACS) 

0.150 

1.2 

0.2 

2012 Dec 01 

MWA 

Hurlev-Walker et al. (2014) (MWACS) 

0.180 

1.9 

0.1 

2012 Dec 01 

MWA 

Hurley-Walker et al. (2014) (MWACS) 

0.235 

3.6 

0.5 

2013 Jul 04 

GMRT 

GMRT Cal 

0.408 

6.6 

0.2 

1969-1974 

MOST 

Large et al. (T981) (MRC) 

0.468 

7.0 

0.7 

1965-1969 

Parkes Interferometer 

Wills (1975J 

0.580 

6.6 

0.4 

1965-1969 

Parkes 

Wills (1975) 

0.610 

6.7 

1.1 

2013 Jul 04 

GMRT 

GMRT Cal 

0.635 

7.3 

0.2 

1965-1969 

Parkes 

Wills (1975) 

0.843 

6.5 

0.3 

1987 Nov 07 

MOST 

Campbell-Wilson & Hunstead (1994) 

0.843 

6.4 

0.2 

1998 Apr 05 

MOST 

Mauch et al. (2003) (SUMSS) 

0.960 

6.1 

0.2 

1965-1969 

Parkes 

Wills (1975J 

1.357 

4.6 

0.1 

2012 Jun 15 

ATCA 

This work 

1.384 

4.4 

0.4 

2008 Jun 19 

ATCA 

Randall et al. (2011) 

1.410 

4.70 

0.07 

1965-1969 

Parkes 

Wills (1975) 

1.687 

3.83 

0.09 

2012 Jun 15 

ATCA 

This work 

1.904 

3.43 

0.08 

2012 Jun 15 

ATCA 

This work 

2.106 

3.14 

0.08 

2012 Jun 15 

ATCA 

This work 

2.307 

2.87 

0.07 

2012 Jun 15 

ATCA 

This work 

2.496 

2.7 

0.3 

2008 Jun 19 

ATCA 

Randall et al. (2011) 

2.513 

2.65 

0.06 

2012 Jun 15 

ATCA 

This work 

2.700 

2.49 

0.03 

1965-1969 

Parkes 

Wills (1975) 

2.718 

2.45 

0.06 

2012 Jun 15 

ATCA 

This work 

2.921 

2.28 

0.05 

2012 Jun 15 

ATCA 

This work 

4.680 

1.37 

0.03 

2010 Jun 10 

ATCA 

This work 

4.800 

1.35 

0.07 

2004 Nov 11 

ATCA 

Murphy et al. (2010) (AT20G) 

4.800 

1.25 

0.02 

1993 Sep 24 

ATCA 

McConnell et al. (2012) (ATPMN) 

4.800 

1.3 

0.1 

2008 Jun 19 

ATCA 

Randall et al. (2011) 

4.850 

1.39 

0.07 

1990 Jun 01 

Parkes 

Wright etal. (1994) (PMN) 

4.926 

1.28 

0.03 

2010 Jun 10 

ATCA 

This work 

5.000 

1.1 

0.2 

1965-1969 

Parkes 

Wills (1975) 

5.009 

1.35 

0.08 

1965-1969 

Parkes 

Wills (1975') 

5.145 

1.22 

0.03 

2010 Jun 10 

ATCA 

This work 

5.380 

1.75 

0.03 

2010 Jun 10 

ATCA 

This work 

5.615 

1.09 

0.03 

2010 Jun 10 

ATCA 

This work 

5.822 

1.04 

0.03 

2010 Jun 10 

ATCA 

This work 

6.062 

0.98 

0.02 

2010 Jun 10 

ATCA 

This work 

6.264 

0.94 

0.02 

2010 Jun 10 

ATCA 

This work 

8.215 

0.65 

0.02 

2010 Jun 10 

ATCA 

This work 

8.478 

0.62 

0.01 

2010 Jun 10 

ATCA 

This work 

8.640 

0.57 

0.06 

2008 Jun 19 

ATCA 

Randall etal. (2011) 

8.640 

0.55 

0.03 

1993 Sep 24 

ATCA 

McConnell et al. (2012) (ATPMN) 

8.640 

0.61 

0.03 

2004 Nov 22 

ATCA 

Murphy et al. (2010) (AT20G) 

8.690 

0.60 

0.02 

2010 Jun 10 

ATCA 

This work 

8.936 

0.58 

0.02 

2010 Jun 10 

ATCA 

This work 

9.156 

0.56 

0.01 

2010 Jun 10 

ATCA 

This work 

9.360 

0.55 

0.01 

2010 Jun 10 

ATCA 

This work 

9.575 

0.54 

0.01 

2010 Jun 10 

ATCA 

This work 

9.776 

0.52 

0.01 

2010 Jun 10 

ATCA 

This work 

18.500 

0.18 

0.01 

2002 Mar 21 

ATCA 

Ricci et al. (2004) 

19.904 

0.152 

0.008 

2004 Oct 22 

ATCA 

Murphy et al. (2010) (AT20G) 

22.000 

0.13 

0.01 

2001 Jan 02 

ATCA 

Ricci et al. (2006) 
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(120141) only reports a flux density at 0.180 GHz based on a 
spectral fit to the flux density at the three frequencies. Us¬ 
ing all three flux density measurements is justified due to the 
brightness of PKS B0008-421 relative to other sources in the 
catalog. Confusion also has a negligible impact on the flux 
density measurements at these frequencies since PKS B0008- 
421 is at least four times brighter than any other source in the 
synthesized primary beam. 

2.3. Archival Data 
2.3.1. ATCA 

PKS B0008-421 has been observed by the ATCA as part 
of both survey and t argete d programs. For targeted observa¬ 
tions, (RandajreFal] (j^Hl) observed PKS B0008-421 as part 
of a campaign to construct an unbiased southern hemisphere 
catalog of CSS and GPS sources. They conducted observa¬ 
tions at 1.4, 2.3, 4.8 and 8.6 GHz in the 1.5B array using 
the now superseded 128 MHz array backend. The source was 
not included as part of their final catalog since the observed 
frequencies were not low enough to constrain the turnover. 
PKS B0008-421 was also observed by iRicci et al.i (120041) as 
part of a program focusing on the polarization properties of 
bright radio sources at 18.5 GHz. Since the 15 mm receivers 
were in the commissioning phase at the time of observation, 
only three antennas in a compact confi gurati on were used with 
128 MHz bandwidth. Additionally, Ricci et all (120041) had 
poor atmospheric phase stability during the observations so 
they used non-imaging model-fitting techniques to calibrate 
and derive source flux densities in all the Stokes paramete rs 
after a ccounting for confus ion. The sources of IRicci et ail 
(120041) were followed up bv IRicci et all (120061) at 22 GHz. 

During the _blind Australia Telescope 20 GHz Survey 
(AT20G: lMurphv et alMTol) . PKS B0008-421 was observed 
with the ATCA using two 128 MHz bands centered on 18.752 
and 21.056 GHz, which were then combined into a single 
band centered on 19.9 GHz. The source was followed up at 
4.8 and 8.6 GHz. Upper limits were placed on the polarized 
flux and fractional polarization at these frequencies. All ob¬ 
servations were conducted in a hybrid array. PKS B0008- 
421 was also observed as part th e Australia Telescope - 
Parkes-MIT-NRAO survey (ATPMN: lMcConnell et al.i2012l) . 
which was an AT CA survey of Parkes-MIT-NRAO (PMN; 
IWright et al.ll 1994h sources at 4.8 and 8.64GHz, with the ar¬ 
ray in 6C and using a 128 MHz bandpass. PKS B0008-421 
was spatially unresolved for all observations conducted by the 
ATCA. Furthermore, PKS B0008-421 remained unresolved 
when the AT20G data_were reprocessed utilizing the study of 
6 km baselines (IChhetri et al.ll2013l) . All ATCA flux densities 
are listed in Table |T] 

2.3.2. Parkes 64 m Radio Telescope 

The Parkes 64 m radio telescope observed PKS B0008-421 
between 1965 and 1969 at 0.580, 0.635, 0.960, 1.410, 2.700, 
5.000 and 5.009 GHz (I Willsl 1T9751 and references therein). 
At 0.580, 0.635 and 0.960GHz uncertainty due to confusion 
can be large since the Parkes 64 m radio telescope is a single 
dish, so scans were made over a large range of right ascen¬ 
sion and declination to correct for effects of any strong con¬ 
fusing sources on the base level on either side of the scan. 
PKS B0008-421 is located in a field that is isolated from 
sources brighter than 1 Jy at 0.635 GHz. PKS B0008-421 was 
also observed by the Parkes interferometer, which consisted 
of the 64 m radio telescope and a 18 m radio telescope, at 


0.468 GHz. This had the advantage of allowing the effects of 
confusing sources on source structure to be separated and re¬ 
moved. These observations are discussed at greater length by 
IWillsI (119751) . Note that the fl ux densities h ave b een adjusted 
from t he flux density scale of lWillsI (119731) to the lBaars et alJ 
(119771) flux density scale via appropriate frequency dependent 
multiplicative factors. Additionally, PKS B0008 -421 was ob- 
served as part of the PMN survey at 4.85 GHz (IWright et alJ 
119941) . The Parkes 64 m radio telescope flux density measure¬ 
ments are also listed in Table [T] 

2.3.3. Molonglo Observatory Synthesis Telescope 

PKS B0008-421 was observed by the MOST during 
the Sydney University Molonglo Sky Survey (SUMSS; 
iBock et alJTl999t iMauch et all 12001 at 0.843 GHz and the 
Molonglo Reference Catalogue (MRC; lLarge et afl 119811) 
at 0.408 GHz. Additionally, PKS B0008-421 was fre¬ 
quently observed by the MOST at 0.843 GHz as part 
of a_ cal ibrat or obse r vatio n ca mpaig n fro m 1 98 7 to 1 996 
(ICampbell-Wilson & Hunsteadl 11994c iGaensler & Hunsteadl 
120001) . The flux density measurements are p rese nted in Ta¬ 
ble HI See the discussion on variability in § 14.11 for further 
details about these observations. 

2.3.4. Giant Metrewave Radio Telescope 

The GMRT observed PKS B0008-421 at 0.235 and 
0.610GHz as part of its observing program to identify suit¬ 
able low radio frequency calibrator^]. PKS B0008-421 was 
unresolved at both frequencies. Uncertainties on the flux den¬ 
sity measurements were calculated from the noise in the im¬ 
ages and the errors in the gain solution of PKS B1934-638. 
The flux densities are presented in Table |T] 

2.3.5. Very Long Baseline Interferometry Observations 

VLBI observations of GPS sources are integral to under¬ 
standing their nature because it is only with this high level 
of resolution that the source structure can be characterized. 
PKS B0008-421 was part of a VLBI monitoring campaign of 
southern hemisphere GPS sources at 2.3 and 8.4 GHz from 
1990 to 1993 using a combination of antennas at s eve n differ¬ 
ent l ocati ons in Aus tralia and one in South Africa (iKingl 1994k 
LTauncev et alJl2003l) . Based on model fitting to the visibili¬ 
ties and closure phases, PKS B0008-421 was well described 
by a model comprising two aligned linear components (IKingl 
119941) . At 2.3 GHz, the separation of these two components 
was observed to be ~ 120 milliarcseconds, which at z = 1.12 
corresponds to a physical separation of ~ 1000 pc; calculated 
assuming a spatially flat ACDM cosmology with matter den¬ 
sity S7 m = 0.286, vacuum energy density Q\ = 0 .714, and 
Hubble constant H 0 = 69.6kms _1 Mpc~ x dWrightl 120061) . 
These structural properties are als o ev ident in VLBA mea¬ 
surements of the source at 8.6 GHz dFev et al.l2004l) . although 
detailed model fitting was not conducted for these observa¬ 
tions. This classical double structure and small separation 
makes PKS B0008-421 consistent with the morphological 
CSO classification of GPS sources. 

3. BAYESIAN INFERENCE MODEL FITTING ROUTINE 

We have implemented a Bayesian inference model fitting 
routine to assess the extent to which different absorption mod- 

1 Details on the data reduction process can be found at 
http://gmrt.ncra.tifr.res.in/gmrt_hpageAJsers/Help/CAL/Cal-List.html 
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Figure 1. Spectral energy distribution forPKS B0008-421, as described by the data points in Table[T] All points plotted are in the observed frame of PKS B0008- 
421. 


els can describe the spectrum of PKS B0008-421, and to cal¬ 
culate the respective model parameters and associated uncer¬ 
tainties. Suppose that the observed flux densities at various 
frequencies are represented by the data vector D and that we 
want to estimate the values of the parameters 0 in some un¬ 
derlying model of the data M. One possible way to approach 
this problem is to utilize Bayes’ theorem, 


Pr(0 | D, M) 


Pr(D | 0,M)Pr(0 \ M) 
Pr(£> | M) 


( 1 ) 


where Pr(0 | D. M) is the posterior distribution of the model 
parameters, representing the updated belief of the model pa¬ 
rameters given the data, Pr(D | 0, M ) = £(0) is the likeli¬ 
hood function, which is the probability of observing the data 
provided given some model parameters, Pr(0 | M) = 11(0) 
is the prior information associated with the model, and the 
Bayesian evidence Pr(£) | M) = Z is a measure of how well 
the model predicts the data observed. 

The flux density measurements of PKS B0008-421 in Table 
IH are independent since the measurements are conducted by 
different instruments or the receivers had a response that en¬ 
sures independence between channels. We make the justified 
assumption that the uncertainty er„ on each measurement is 
normally distributed. 

Since the flux density measurements are independent and 
the uncertainty on each measurement is normally distributed, 
the joint log likelihood function we used for this analysis was: 


ln£(0) 




(Ai - A(^n )) 2 


+ ln(27rcr^) 


, ( 2 ) 


where D n is the flux density observed at frequency u n and 


S v (v n ) is the expected flux density at v n from the model given 
parameters 0. This likelihood function was convolved with 
uniform priors for each model parameter with sensible ranges 
placed by known physical constraints. For example, the peak 
frequency or normalization of the flux density can not be neg¬ 
ative. 

The most efficient way to calculate the posterior probability 
density functions for the parameters of different models is via 
Markov chain Monte Carlo (MCMC) methods, which are de¬ 
signed to discretely sample posterior probability density func¬ 
tions of model parameters such that the likelihood function 
is maxi mized. We implemented an affine-invariant ensemble 
sampler dGoodma n & Wea rel 120101). via the Python package 
emcee (iForeman-Mackev et al.ll2013h . which means the algo¬ 
rithm utilizes an ensemble of “walkers” and is relatively in¬ 
sensitive to covariance between parameters since it performs 
equally well under all linear transformations. We utilized the 
simplex algorithm to direct the walkers to the maximum of 
the likelihood function dNelder & MeadlU965l) . 

A Bayesian inference approach to modeling radio spec¬ 
tra has several advantages over traditional probabilistic data 
analysis procedures, such as least-squares and maximum- 
likelihood methods. This includes the ability to marginal¬ 
ize over nuisance parameters, generate complete probability 
distributions of model parameters and to incorporate prior 
knowledge so known physical constraints can be naturally 
placed on the system. Most importantly, since we are fit¬ 
ting non-linear models, Bayesian inference facilitates a more 
objective model comparison than traditional reduced x‘ 2 ~ 
tests via comparison of evidence values, which appropri¬ 
ately penalize any additional degrees of freedo m by i n tegrat- 
ing over all param eter space in each model (iJeffrevsl 1 196 lL 
lAndrae et al.li2010t) . 
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Bayesian evidence is required to normalize the posterior 
over the prior volume so that: 

z = J J ■■■ j c(e)u(6)de, ( 3 ) 

where the dimensionality of the integration is set by the num¬ 
ber of parameters in the model. The evidence can be thought 
of as the numerical statement of Occam’s razor, since it repre¬ 
sents the average of the likelihood over the prior volume. This 
means that models with high likelihood values throughout pa¬ 
rameter space are favored over models with low likelihood 
regions. Unless a complex model provides a significantly bet¬ 
ter fit to the data than a model with a smaller parameter space, 
the evidence value will be larger for the model with a smaller 
number of parameters. We calculated the evidence for the dif¬ 
ferent models using the algorithm MultiNest dFeroz et al.l 
120131) . which is an implementation of nested sampling. The 
algorithm is initialized by uniformly sampling the prior space 
and then contracting the distribution of samples around high 
likelihood regions by discarding the sample with the least 
likelihood. A random sample is duplicated to keep the num¬ 
ber of samples constant. This process is repeated until the re¬ 
placement of samples is optimized after finding the region of 
maximum likelihood. The strength of nested sampling is that 
it calculates the mean posterior probability and the evidence. 
We ensured that the parameter estimates from MultiNest 
were within uncertainty of those outputted by emcee, such 
that the calculated evidence was indicative of the best model 
fit. 

Assuming a priori that two models are equally likely to 
describe the data, model selection can be performed solely 
based on the ratios of the respective evidences. In this paper 
we perform this in log space such that Aln(Z) = ln(^ 2 ) — 
lnf-Zj), Using an updated version of the Je fferys scale (e.g. 
IKass & RaftervHT995l:fScaJt'e & Healdl201 21) . A In (Z) > 3 is 
strong evidence that M 2 is favored over Mi, 1 < Ain (Z) < 
3 is moderate evidence in support of M 2 over M\, while 
0 < Ain (Z) < 1 is inconclusive. For the sake of compar¬ 
ison to literature we have also included reduced x' 2 values for 
model selection assessment. 

Note that each flux density measurement will have a sys¬ 
tematic uncertainty associated with the absolute flux density 
scale, which could cause the ensemble of measurements to de¬ 
viate from a Gaussian distribution. To test whether this would 
have an impact on the fitting proc ess w e introduced hyper¬ 
parameters to the fits dHobson et al.H2002[) . Hyperparameters 
are similar to nuisance parameters in that they are marginal¬ 
ized over when calculating the posterior distribution, but hy¬ 
perparameters differ in that they are not present in the model 
beforehand. They provide a method to quantify inaccurate un¬ 
certainties and syste matic errors via weighting of the data sets 
dHobson et al.l2002l) . We find no statistical evidence that a hy¬ 
perparameter is justified to model a possible deviation from a 
normal distribution for the ensemble of measurements, with 
Ain (Z) < 1 when comparing fits with and without a hy¬ 
perparameter. The reason the ensemble distribution closely 
adheres to a Gaussian distribution is likely due to the fact that 
PKS B1934-638 was used to set the flux density scale for the 
majority of the observations. 

4. RESULTS 

The different absorption models for radio galaxies can be 
separated into two broad categories based on the underlying 


absorption mechanism: FFA and SSA. For FFA, the models 
vary depending on whether the ionized screen is internal or 
external to the radio emitting plasma, and on the topology 
of the screens. Additionally, the spectrum can also exhibit 
breaks due to the ageing of the electron population. In this 
section we first justify the use of multi-epoch data to fit the 
spectra of PKS B0008-421, and then outline the different fea¬ 
tures of the absorption models and fitting statistics. 

4.1. Variability 

When considering the use of multi-epoch data to de¬ 
scribe the spectrum of PKS B0008-421, it is important to 
determine whether PKS B0008-421 displays any tempo¬ 
ral variability that could result in deviations to its intrin¬ 
sic spectrum. PKS B0008-421 has been part of several ra¬ 
dio monitoring campaigns, since it is used as a primary 
or secondary cal ibrato r for many so uthern hemisphere ra- 
dio tel esco p es. ICampbell-Wilson & Hunsteadl (119941) and 
iGaensler & Hunsteadl OOPOT both studied PKS B0008-421 
with the MOST to ascertain whether it displayed any evi¬ 
dence of intrinsic variability or scintillation at 0.843 GHz. 
PKS B0008-421 was observed irregularly over 22,000 times 
in a 2-minute scan mode between 1987 and 1996, with sep¬ 
aration between the observations ranging from several hours 
to days. In order to test for variability. IGaensler & Hunsteadl 
(120001) calculated the x' 2 probability that the flux density has 
remained constant from the light curve of PKS B0008-421, 
binned in 30 day intervals. This bin size was selected to re¬ 
duce the effects of unrecognized systematic errors in the flux 
density determination and the presence of confusing sources 
in the field. PKS B0008-421 was found to be non- variab le 
on this time scale following the procedure of lKesteven et all 
(119761) . which compares the probability of whether the disper¬ 
sion in the weighted mean flux density comes from intrinsic 
variability or distributions of measurement errors, with a p- 
value significantly greater than 0.01. 

PKS B0008-421 was also part of a monitoring campaign by 
the Tasmanian Mount Pleasant 26 m antenna between 1990 
and 1993. PKS B0008-421 was regularly observed over a 30 
month period 52 and 34 times at 2.3 and 8.4 GHz, respectively 
(iJauncev et al .1120031) . The average root-mean-square errors 
of the observations were ~ 0.1 Jy. ..Using the same statisti¬ 
cal technique applied by IGaensler & Hunsteadl (120001) . I Kind 
(119941) demonstrated that PKS B0008-421 was not variable at 
2.3 GHz and 8.4GHz with p-values of 0.02 and 0.95, respec¬ 
tively. 

Therefore, it is evident that PKS B0008-421 is not variable, 
within the uncertainties of these datasets, over timescales 
from hours to decades at several frequencies. This lack of 
variability is consistent with the nature of most GPS sources 
that have morphologies c onsistent with CSOs (lO’Dealll998l : 
iFassnacht & Tavlorl2001l) . and justifies the use of multi-epoch 
data to model the spectrum. 

4.2. Free-Free Absorption 

FFA involves the attenuation of emission radiation by an ex¬ 
ternal or internal ionized screen relative to the emitting elec¬ 
trons. The morphology of the screen could be either homo¬ 
geneous or inhomogeneous. For electron temperature T e in K 
and free electron density n e in cm" 3 , the optical depth due to 
FFA at frequency v in GHz is approximated by 

t v w 8.24 x 10 -2 v -2 - 1 T“ 1 - 3B J n\ d l, 


(4) 
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where l is the distance through the source in pc 
(iMezger & Hendersonl 119671) . Assuming that a homoge¬ 
neous ionized screen surrounds the plasma producing the non- 
thermal power-law spectrum, the resulting free-free absorbed 
spectrum is described by 

S v = av°‘e~ T ’', (5) 

where a characterizes the amplitude of the intrinsic syn¬ 
chrotron spectrum, and a is the spectral index of the syn¬ 
chrotron spectrum. It is convenient to parameterize the optical 
depth as tv = (v/v P ) -21 , where v p is the frequency at which 
the optical depth is unity. The best model fit for the spectrum 
is shown in Figure 0a) and the resulting parameter values, 
with associated 1-er uncertainties, are presented in Table 0 

As is clear from the fitting statistics, a single homogeneous 
absorbing screen cannot accurately describe the spectrum 
of PKS B0008-421: the gradient below the turnover is too 
steep and it overpredicts the flux density at high frequencies. 
One p ossible alternative, as suggested by iTingav & de Kooll 
<0X1), is to fit two separate homogeneous free-free absorb¬ 
ing screens in front of two separate non-thermal distributions 
of electrons, since the source has been resolved into two com¬ 
ponents. The justification for this method is that for all obser¬ 
vations with lower resolution than VLBI, the flux density we 
measure represents the sum of the environments of both these 
components. Therefore, the double homogeneous free-free 
model fit to the spectrum has the form: 

S„=Y, ( 6 ) 

{= 1,2 

The resulting fit is also presented in Figure 0a) and Table 0 
The double homogeneous FFA model is a statistical improve¬ 
ment over the single homogeneous model, with A In (Z) > 6 . 
However, despite doubling the number of parameters in the 
model, the double homogeneous FFA model still falls too 
rapidly at low radio frequencies and over predicts the flux den¬ 
sity at high frequencies. 

It is possible that the absorbing ionized medium is mixed 
in with the relativistic electrons that are producing the non- 
thermal spectrum. In that case the spectrum is characterized 
as 



This fit is presented in Figure0b) and in Table 0 It is evident 
that the internal FFA cannot accurately describe the spectrum 
because of its shallow slope below the turnover. This is to 
be expected since the gradient below the turnover for such a 
model should be a + 2.1. 

An alternative FFA model, proposed by iBicknell et al.l 
(IT997h . treats the FFA screen as inhomogeneous and external 
to the non-thermal electrons in the lobes of the source. In this 
model, the jets from the AGN produce a bow shock that pho- 
toionizes the interstellar medium as it propagates, such that 
if the density of the gas surrounding the radio lobe is suffi¬ 
ciently high, it will ca use FFA . Hence, the clouds have a range 
of optical depths that IBicknell et al.l (119971) assume follow a 
power-law distribution r/ parameterized by the index p, such 
that rf oc f (n 2 T~ 1 - 35 ) p dl. Note that we require p > — 1 
otherwise the model reduces to the homogeneous limit, and 
rj varies from 0 to a maximum value rj 0 = z/“ 21 . Assuming 
that the length scale of the shock and the inhomogeneities in 


the interstellar medium are both much less than the size of the 
lobes, the model can be represented by 


S„ 


a(p+ 1)7 




2.1(p+l)+a 

, ( 8 ) 


where 7 is the lower incomplete gamma function of order 
p + 1. The model is plotted in Figure 0c) and the fit pa¬ 
rameters are represented in Table 0 This is the best fitting 
FFA model with A In (Z) > 3 when compared to the double 
homogeneous model. Comparatively to the double homoge¬ 
neous FFA model, the inhomogeneous FFA model provides 
a better fit to the slope below the turnover but continues to 
predict additional flux at high frequencies. 


4.3. Synchrotron Self-Absorption 

SSA has been the other dominant theory used to explain 
the inverted spectra of GPS sources. In this model, the 
turnover occurs because the source cannot have a brightness 
temperature that exceeds the pl asma temperature of the non- 
thermal electrons dKellermannll 19666 . Alternatively, SSA can 
be thought to occur at the frequency at which the relativistic 
electrons and emitted synchrotron photons have a high chance 
of scattering. The absorption cross-section for a synchrotron 
electron and a low-energy photon is greater at larger wave¬ 
lengths. Hence, for a source of a given size, very-long wave¬ 
length emission is only visible from a very thin shell at the 
surface of the source. As the observing frequency increases, 
photons emerge from regions of the source that are progres¬ 
sively deeper in the source such that the total flux density in¬ 
creases until the optically thin regime is reached. Parameter¬ 
izing the spectrum in terms o f the power-law index, 0, of the 
electron energy distribution (ITingav & de Koof 12003 ). such 
that a = — (ft — l)/ 2 , and assuming the synchrotron emitting 
region is homogeneous, the spectrum is modeled by: 


where. 


S v — (X 


fer i>/2 


1 — e 


r = 




(9a) 


(9b) 


and where is p corresponds to the frequency at which the 
source becomes optically thick. This is the frequency at which 
the mean free path of electron-photon scatterings is approxi¬ 
mately the size of the source. Applying the same logic as 
in §[40 that the observed flux density from the source is the 
sum of two different lobes with different electron energy dis¬ 
tributions, the spectrum could also be described by two SSA 
components such that 


S v 


H a i 

i=l,2 


(v_ \- ( ^- 1)/2 


1 - e~ Ti 


Ti 


where. 


Ti = 


V > i l 


-(/3i+4)/2 


( 10 a) 


( 10 b) 


The fits for the single and double SSA models are presented in 
Table0and Figure0d). The double SSA model is statistically 

























Table 2 

Best fit parameters and their associated uncertainties for the different absorption models. The uncertainties presented are 1 -<t. Note the reported values are for the observed frame of PKS B0008-421. Free-free 
absorption and synchrotron self-absorption are abbreviated to FFA and SSA, respectively. The parameters listed in the table are: the normalization parameters of the intrinsic synchrotron spectrum a;, the spectral 
indices of the synchrotron spectrum ctj, the power-law indices of the electron energy distribution /3j, the frequencies at which the optical depth of the absorption model is unity v p ,i, the high-frequency 

exponential cutoff v\, r , the reduced x 2 -' , alue of the model fit x 2 ed . anc * 1°S °f the Bayesian evidence of the model fit ln(Z). 
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better fit than the single SSA model, with Ain (Z) > 100. 
However, as evident in Figure Efd), both models accurately 
model the slope below the turnover but underestimate the flux 
density at low frequencies and overestimate it at high frequen¬ 
cies. 

We do not outline an inhomogeneous model of SSA since 
the slope below the turnover in the spectrum of PKS B0008- 
421 is close to the theoretical prediction of the homogeneous 
SSA model. An attempt to fit a more general inhomogeneous 
SSA model finds that the steep spectral slope forces the model 
fit to be degenerate with the homogeneous SSA model. 

4.4. Spectral breaks and cut-offs 

As the non-thermal electrons age in the jets and lobes of a 
radio source, synchrotron and inverse-Compton losses pref¬ 
erentially deplete high-energy electron populations such that 
the spectrum deviates from a single power-law. There are sev¬ 
eral models that predict how the shape of the spectr um s hould 
evolve. In the continuous injection model of iKardashevi 
(119621) a stream of electrons is constantly injected into a vol¬ 
ume permeated with a constant magnetic field. This model 
predicts that the optically thin part of the spectrum should 
steepen above a break frequency, v' hl , from a to a — 0.5. Once 
this characteristic break frequency is known, the sourc e’s ag e 
and magnetic field strength can be estimated (iMurgia et alJ 
119991) . We find that augmenting the models in the previ¬ 
ous section with the continuous injection model does not 
improve the fit to the data because of its discontinuous na¬ 
ture. The abrupt change in spectral index is inconsistent with 
the smooth, continuous decline evident in the three bands of 
CABB data in Figure [Q and the fact the decrease in slope is 
not steep enough to fit this decline. 

For GPS sources the assumptions of the continuous injec¬ 
tion model may not hold, since there is evidence of intermit- 
tent and recurrent activity i n the nuclei of such so urces (e.g. 
lOwsianik & Conwavlll998l:lBrocksopt> et al.ll2007l) . Approxi¬ 
mately ten percent of GPS sources display several large scale 
radio structure s that are usually explain ed as relics of previous 
active phases (IStanghellini et al.|j2005 ). An alternative model 
is to assume that the injection of fresh particles has ceased 
and that the radio source is already in the relic phase. In ad¬ 
dition to the transition to the steeper slope a — 0.5 after u^ T , 
this model predicts the formation of a high-frequency expo¬ 
nential cutoff after a second break frequency, such that 
the absorption models have an additional multiplicative factor 
e -^Abr (|j a ffg & Perolall 1 97.1t iKomissarov & Gubanovll 19941 : 

lMurgiall2003l) . This second break frequency is related to the 
first according to 


Ubr = z-C > C 11 ) 

where t s is the age of the source and t D g is the time since the 
injection of fresh electrons has ceased. 

As mentioned above, a steepening of the optically thin slope 
to a — 0.5 is not statistically evident for PKS B0008-421. 
However, it is possible that the first break frequency iy' hr has 
moved to a lower frequency below the spectral turnover. This 
is consistent with the fact that the absorption models require 
an extraordinarily steep injection spectral index of a ~ —1.2 
to fit the high-frequency data of PKS B0008-421. If the in¬ 
jection has ceased and the first spectral break frequency has 
moved to the optically thick part of the spectrum, a more rea¬ 
sonable injection spectral index of a ~ —0.7 can fit the opti- 
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Figure 2. Different absorption model fits to the spectral energy distribution of PKS B0008-421, using the parameter values reported in Table [2] The colored 
regions associated with the models represent the l-cr uncertainty on the model fit at the respective frequency. The y-values for the model fits are displayed below 
the spectral energy distributions. Note that free-free absorption and synchrotron self-absorption in these figures are abbreviated to FFA and SSA, respectively, (a) 
Single and double homogeneous FFA models are shown in brown and red, respectively, (b) Internal FFA model is represented in gray, (c) Inhomogeneous FFA 
model is shown in green, (d) Single and double SSA model fits are presented in purple and blue, respectively, (e) Inhomogeneous FFA model with an exponential 
spectral break is shown in green, (f) Single and double SSA models, with exponential spectral breaks, are presented in purple and blue, respectively. All fits are 
conducted in the observed frame. 





















































































































































Spectral Modeling of GPS Source PKS B0008-421 


11 


cally thin part of the spectrum. As is clear in Table [2] the addi¬ 
tion of an exponential attenuation at high frequencies provides 
a very significant statistical improvement in the overall fit to 
the spectrum, with Ain (Z) > 100 for corresponding mod¬ 
els compared with and without the exponential break. Note 
that while all the models are improved by the exponential 
break, only the three best fitting models are presented in Fig- 
ures[2}e) and[2jf) as the addition of the exponential break still 
does not allow the spectral slope below the turnover to be well 
fit for homogeneous free-free and internal FFA models. Af¬ 
ter the addition of the high-frequency exponential break, the 
two best fitting models are the double SSA and the inhomo¬ 
geneous FFA. The difference in Bayesian evidence between 
these two models, Ain (Z) = 1.7, moderately favors the in¬ 
homogeneous FFA model with an exponential break over the 
double SSA model with an exponential break. Such expo¬ 
nential spectral breaks have been obs erved in other sources, 
for instance PKS B1518+047 (lOrienti et al.tl2010t) . Signifi¬ 
cantly, not only does this addition allow the models to follow 
the curvature evident in the CABB bands, it also ensures that 
the double SSA and inhomogeneous FFA models do not un¬ 
derestimate the flux density at low frequencies. This indicates 
that it is necessary to have good high-frequency coverage to 
correctly model the flux below the turnover in the spectrum. 

Note that the high-frequency steepening could be due to 
source structure gradually being resolved out when observing 
from low to high frequencies or by the presence of older ex¬ 
tended structure. As discussed in §[2] the steepening cannot 
be due to resolution affects since PKS B0008-421 was unre¬ 
solved at all frequencies, even when 6 km baselines are used 
at 20 GHz (IChhetri et al.ll2013l) . All imaging at the frequen¬ 
cies of the ATCA observations, and at VLBI resolutions, also 
did not provide any evidence for extended structure. There¬ 
fore, the break evident in the spectrum is highly likely to be a 
physical feature of PKS B0008-421. 

5. DISCUSSION 

PKS B0008-421 represents an extreme outlier in the pop¬ 
ulation of GPS sources. It has the smallest spectral width 
(defined as the fullwidth at half maximum), ~ 0.5 decade of 
frequency, and steepest low radio frequency slope, a ~ +2.4, 
of any known GPS source. The average spectral width and 
spectral index below the turnover for GPS sources are ~ 1.2 
decades_of frequency (lO’Dea et all fl99ll) and a ~+0.7 
(Ide Vries et all Il997l) . respectively. The outlier nature of 
PKS B0008-421 is emphasized in Figure [3] which shows 
the distribution of the spectral indices of MWACS sources 
that have a AT20G counterpart. Such a population is domi¬ 
nated by bright, compact radio sources, likely including many 
GPS sources. The spectral slope of PKS B0008-421 in the 
MWACS band is over five standard deviations away from the 
median of the distribution. 

Since PKS B0008-421 has an observed peak frequency be¬ 
low a gigahertz and a double component morphology that has 
a physical separation of ~ 1000 pc, it has several similari¬ 
ties to a CSS source. However, when the relationships be¬ 
tween the observed spectral peak, radio power, angular size 
and p eak bri ghtness are compa red (e.g. Fig. 1. ISnellen et all 
12000: Fig. l. lAn &Baanll2012l) . all of PKS B0008-421 spec¬ 
tral properties are consistent with the GPS/CSO classifica¬ 
tion. Additionally, the physical separation between the double 
components is too large to for PKS B0008 -421 t o be associ¬ 
ated with the HFP population (lOrienti et al .1120061) . 

It is noteworthy that the slope of the spectrum below the 



“MWACS 

Figure 3. Distribution of the MWACS spectral index for the 706 sources 
detected in both MWACS and AT20G. Such a sample is dominated by com¬ 
pact sources, a in this plot represents the spectral index of the source in the 
MWACS observations only. The location of PKS B0008-421 is highlighted 
in orange and by an arrow. This location is over five standard deviations away 
from the median of the distribution. 

turnover of PKS B0008-421 is close to the SSA theoreti¬ 
cal limit of a =+2.5. This is because the common ex¬ 
planation for observing shallower spectral gradients below 
the turnove r is due to inhomogeneity of the SSA regions in 
the source (IQ’Deal 119981 iTingav & de Kooill20()3h . Such a 
steep low-frequency gradient could be indicative of a homo¬ 
geneous emitting structure that could allow the most compre¬ 
hensive differentiation between the different types of absorp¬ 
tion models. Additionally, the very steep optically thin spec¬ 
tral index of a 1.2 is a surprising characteristic since, if 
the youth scenario is correct, the radio emission from GPS 
sources should have only started several thousand years ago. 
This would imply an extraordinary injection spectral index is 
needed, or a source significantly older than common for GPS 
sources. 

As shown in §|4j we can statistically conclude that inter¬ 
nal FFA, single and double homogeneous FFA, and single 
SSA cannot be responsible for the turnover in the spectrum 
of PKS B0008-421. The two models that best fit the spectrum 
of PKS B0008-421 are the inhomogeneous FFA and double 
SSA models, both with the addition of a high-frequency ex¬ 
ponential break. Both models can accurately predict the flux 
density and slope of the spectrum at high and low radio fre¬ 
quencies. By eye the two fits are identical but the difference 
in Bayesian evidence, A In (Z) = 1.7, moderately favors the 
inhomogeneous FFA model over the double SSA model. 

However, both the inhomogeneous FFA and double SSA 
models require some abnormal physical parameters to gener¬ 
ate their fit, underscoring the extreme nature of PKS B0008- 
421. Using the best fitting parameters from the double SSA 
model with an exponential break, the magnetic field strength 
of the source can be estimated. Provided that the spectral 
turnover of GPS sources is caused by SSA, the turnover fre¬ 
quency axi in GHz, for a homogeneous, self-absorbed ra¬ 
dio source with a_power-law electron energy distribution can 
be expressed (iKellermann & Paulinv-Tothli 1981 ) in terms of 
the magnetic field strength B, in G, the flux density at the 
turnover frequency ,Sj nax , in Jy, and angular size of the source 
0 max , in milliarcseconds, by 

Umax - 8.1 + Z) 1 ' 5 . (12) 
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For PKS B0008-421, this is estimated to be B ~4.1 G. This 
is an order of magnitude larger than the typical magnetic 
field strength reported for GPS sources, which are often in 
the ra nge of 5 to 100 mG (lO’Deall 1 998t lOrienti & Dallacasal 

l2008h . 

In comparison, the inhomogeneous FFA model demands a 
large, positive p value to fit the steep low-frequency spectral 
slope, as shown in Table [2] All previous literature values re¬ 
port values of p that are negative or zero because of the mo re 
common shallow slope below the turnover dBicknell et all 
119971 iTingav & de Kooll l2003h . Since p parameterizes the 
distribution of optical depth of absorbing clouds, the impli¬ 
cation of a large, positive p is that there are substantially more 
clouds of high emission measure than clouds of low emis¬ 
sion measure. In terms of the environment of the source, this 
could be justified if the radio jet is beginning to propagate into 
the environment and breaking up the interstellar medium into 
several large, dense and cool clouds. Therefore, a positive p is 
consistent with the idea that the source is confined to a small 
spatial scale due to high density, low temperature clouds. If 
the jets continue to propagate into the medium, the jet inter¬ 
actions with the surrounding medium will increase the tem- 
perature and decrease the density of the clouds, as evident in 
[Wagner & Bicknel l (1201 If) , such that p will decrease towards 
more common literature values. 

Independent of the underlying absorption model, the ne¬ 
cessity of a high-frequency exponential break to accurately 
fit the spectrum means that injection or accele ratio n of rela¬ 
tivistic particles has ceased in the source (lOrienti et al.ll201flh . 
The lack of variability also suggests that PKS B0008-421 is 
no longer being powered. The presence of a high energy ex¬ 
ponential break, but no statistical evidence of the first spectral 
break frequency, indicates that the first break has moved to the 
optically thick part of the spectrum, v' hr < 0.59 GHz. This al¬ 
lows a more typical injection spectral index a ~ —0.7 to fit 
the optically thin part of the spectrum. With the knowledge 
of the spectral break frequencies, we can place a lower limit 
on the age the source t s and calculate the time since the last 
injection of electrons t Q s using 

t s = 5.03 x 10 4 £r 3/ X r (l + z) 1/2 yr, (13) 

and Equation!]!] where B is in mG and the break frequencies 
are in GHz. 

We can estimate the magnetic field strength of PKS B0008- 
421 independently of SSA through equipartition theory, as¬ 
suming that the total energy densities of cosmic rays and 
the magnetic fields are equal. Using the revised formula of 
[Beck & Krause! (120051) . the equipartition magnetic field for an 
injection spectral index of a = —0.7 is B eq ~6.2 mG. This 
magnetic field strength is inconsistent with that derived from 
SSA theory, but consistent with other magnetic field s trengths 
derived for GPS sources (lOrienti & Dallaca~sall2008 f). Using 
the equipartition magnetic field strength, the radiative lifetime 
of PKS B0008-421 is calculated to be t s > 2900 years and, 
using t'br from the inhomogeneous FFA model, the time since 
the last injection of fresh electrons is t 0 s = 550 ±110 years. 
Both of these value s are c onsistent with the youth hypothe¬ 
sis of GPS sources dMurgia et al.lll9991) . Applying the much 
larger magnetic field strength calculated from SSA theory and 
the fitted ;/( jr from the double SSA model, the radiative age 
lower limit is unrestrictive and we find t 0 s ~ 650, 000 years. 
This time period is inconsistent with a source with such a large 
flux density and provides a physical argument in favor of in¬ 


homogeneous FFA over the double SSA model. 

Further evidence that the inhomogeneous FFA model is re¬ 
sponsible for the turnover in the spectrum of PKS B0008-421 
is that the peak frequency occurs at around a gigahertz. If 
the injection of relativistic electrons has finished, strong adia¬ 
batic cooling should see a shift of the peak frequency quickly 
out of the gigahertz regime. Assuming the magnetic field is 
frozen into the plasma, the spectral peak would decrease from 
v n n at time tn to v n i at t ime t\ with zz Pi i = z'p.o^oAi) 4 
(lOrienti & Dallaca~sa 12008 ). This would suggest the detec¬ 
tion of a GPS source that has ceased actively fueling radio 
jets would be rare because the peak would quickly shift out of 
common observable radio frequencies. However, a possible 
way that PKS B0008-421 has maintained its peak frequency 
close to a gigahertz, despite the cessation of fresh electrons 
around 500 years ago, is through the presence of a dense 
ambient medium. This dense medium would limit the adi¬ 
abatic losses and maintain the observed peak brightness and 
frequency. This again favors the inhomogeneous FFA model 
to explain the spectrum of PKS B0008-421 since the large, 
positive p parameter needed for the fit also implies a dense 
surrounding medium. A test of this hypothesis would be to 
target PKS B0008-421 for Hi absorption studies and to de¬ 
termine whether the associated absorption towards the source 
has a column density > 10 2O cm -2 , since this model sug¬ 
gests it is pr obable much of the medium surrounding the je ts 
is unionized (IVermeulen et al.ll2003HTengstrand et al.ll2009h . 

It is important to emphasise that both the double SSA and 
inhomogeneous FFA model fits require extreme physical con¬ 
ditions to provide an adequate fit to the spectrum PKS B0008- 
421. In the case of the double SSA, the model requires a 
magnetic field strength that is an order of magnitude larger 
than the average for the population of GPS sources. While for 
the case of FFA, the model requires a distribution of emission 
measure that is dominated by cold and dense clouds to pro¬ 
vide an adequate fit. However, the magnetic field strength re¬ 
quired by SSA is inconsistent with equipartition and with the 
characteristic time scales required by the exponential spectral 
break. Since the physical parameters required by FFA are in¬ 
dependently consistent with the characteristic time scales of 
the exponential spectral break and with the fact that the peak 
of the spectrum is still near a gigahertz, we suggest that the 
extreme physical environment that the inhomogeneous FFA 
model requires is more likely than that required by the double 
SSA model fit. 

The study of PKS B0008-421 also has implications for the 
types of sources that will be discovered by the all-sky surveys 
being performed MWA and LOFAR. Many GPS sources that 
are at z > 1 and have ceased actively fueling their jets would 
not be identifiable at gigahertz wavelengths as classical GPS 
sources, since their spectral peaks would be below 0.3 GHz. 
This means that these two telescopes could reveal high red- 
shift GPS sources and a whole new population of GPS sources 
that have ceased activity and are fading. If a large number of 
GPS sources are identified as having ceased activity, it will be 
possible to conduct a comprehensive analysis as to whether 
such a population can explain the overabundan ce p roblem of 
CSS and GPS sources e ncountered bv lReadhead et al.l (119961) 
and I An & Baanl (120121) . PKS B0008-421 also has implica¬ 
tions for the identification of ultra-steep spectrum sources 
(US S), which are sources defined as having a < —1. While 
iKlamer et al, C 2006I) did not find high-frequency steepening 
for a sample of USS sources, PKS B0008-421 would be 
identified as a USS source if no observations had been con- 
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ducted below 0.6 GHz. This could imply that at least part 
of the USS population consists of dying GPS sources where 
the spectral break has shifted below the low radio frequency 
turnover. Finally, the next generation of blind, large HI ab¬ 
sorption studies, such as the planned HI absorption survey on 
the Australian Square Kilometre Array Pathfinder (ASKAP; 
Ijohnston et al .120(481) . will be vital in constraining the HI den¬ 
sities in GPS sources and could provide additional evidence 
for FFA. 

6. CONCLUSIONS 

We have presented broadband spectral modeling of the ex¬ 
treme GPS source PKS B0008-421, which represents an out¬ 
lier in the current GPS population since it has the steepest 
spectral slope below the turnover, and smallest spectral width, 
of any known GPS source. PKS B0008-421 is part of the CSO 
morphological class, and its lack of variability allowed the use 
of multi-epoch data to describe the spectrum, with spectral 
coverage from 0.118 to 22 GHz. In this analysis we have: 

1. Determined the two best fitting spectral models are dou¬ 
ble SSA and inhomogeneous FFA. We also statistically 
excluded internal FFA, single and double component 
homogeneous FFA, and single component SSA as be¬ 
ing explanations for the turnover in the spectrum of 
PKS B0008-421. 

2. Demonstrated the necessity of a high-frequency expo¬ 
nential break to adequately describe the spectrum of 
the source. This implies that the source has a spectral 
break below the turnover and has likely ceased the in¬ 
jection of fresh particles. With the requirement of a 
high-frequency exponential break, the inhomogeneous 
FFA model is moderately statistically favored over the 
double SSA model. 

3. Highlighted that since activity in the source has ceased, 
constraints can be placed on the time since last injec¬ 
tion. If the physical parameters derived from the dou¬ 
ble SSA fit are used, the time periods that result are 
unrealistic. However, parameters from the inhomoge¬ 
neous FFA model provide a time since the last injec¬ 
tion of fresh electrons t 0 g = 550 ±110 years. The 
existence of a spectral peak in the gigahertz range re¬ 
quires a high density ambient medium to restrict adi¬ 
abatic losses. This conclusion is congruous with the 
dense ambient medium required by the inhomogeneous 
FFA fit. These two lines of evidence lead us to conclude 
that inhomogeneous FFA is responsible for the turnover 
in the spectrum of PKS B0008-421. 

These conclusions imply it is likely that MWA and LOFAR 
will expose a new population of CSS and GPS sources that 
have ceased activity. Such sources would not have been iden¬ 
tified from previous high-frequency surveys since the spectral 
peak would have shifted below a gigahertz. If such a popula¬ 
tion is found to be prominent, it could help explain the over¬ 
abundance of GPS sources in the local Universe relative to 
the number of large galaxies. Additionally, the steep optically 
thin component of PKS B0008-421 suggests it is possible that 
part of the USS population consists of a population of dying 
GPS sources for which the spectral break has shifted below 
the low radio frequency turnover. 
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